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The B+ meson is composed of two heavy quarks of distinct flavor. Measurements of its lifetime and production 
properties have been made based on semileptonic — > J/t/} + l + + X decays using data collected with 
the CDF II detector corresponding to an integrated luminosity of 1 fb — 1 . The B~£ average lifetime cr is 
measured to be 142.5^^'g (stat)±5.5(sj/st) p,m. The measurements of the ratio of the production cross section 

times branching ratio of — > J/rpfj&v relative to B+ — ► J/x/jK+ were done for two P'i (II) thresholds: for 
p T (B) > 4 GeV/c as 0.295 ± 0.040 (stat.)to oil ( s y st -) ± °- 036 (pt spectrum) and for p T (B) > 6 GeV/c as 
0.227 ±0.033 (stat.)±o'.oiT (syst.) ± 0.014 ( PT spectrum). 



1. Introduction 

The meson [1] is composed of an anti-bottom 
quark b and a charm quark c. The presence of two rel- 
atively heavy quarks with different flavors is unique to 
the -B+ system and affects the decay and production 
properties. The theoretically predicted lifetime [2] 
is about a factor of three times smaller than that of 
other B mesons. The expected production cross 
section [3] is about 3 orders of magnitude lower than 
the production cross section of the B + [4] . The first 
observation of the was made using data taken with 
the CDF detector at the Fermilab Tevatron during 
run 1 [5] . Precise mass measurements have been made 
by the CDF Collaboration using fully reconstructed 
£?+ — * J/-07T+ decays, where J/ip decays through 

J/V> -> [6]. 

In this work we report preliminary measurements 
of the B~£ lifetime in the semileptonic decay modes 
J/ipfj, + X and J/ipeX, and the production cross sec- 
tion times branching ratio of the decay mode B+ — ► 
J /tpfj, + v relative to the B + — ► J/ipK + decay. The re- 
sults presented here are based on a data sample with 
an integrated luminosity of 1 fb _1 at -y/s=1.96 TeV 
collected by the CDF II detector. 



2. The B+ lifetime measurement 
concept 

To measure the lifetime of the B+, we construct a 
per event lifetime that is defined using variables mea- 
sured in the transverse plane. If all of the decay prod- 
ucts of the .B+ decay are identified, the lifetime ct is 
the lifetime of the B+ meson in its rest frame mea- 
sured in units of microns of light travel time. It is 



expressed as 



ct 



mL x 
PT 



(1) 



where m is the mass of the B^ , px is the momentum of 
the B+ in the plane transverse to the direction of the 
proton beam, and L xy is the decay length of the B^ 
projected along the transverse momentum. The mass 
of the B+ used in this measurement is m — 6.286 
GeV/c 2 [6]. However, we do not measure all of the 
particles in the semileptonic B+ final state. Instead 
we must define a pseudo lifetime 



ct* 



mL xy (J / j>l + ) 
Pt(J/W + ) 



(2) 



where L xy and px are evaluated using the J/ip + l + 
system. We can obtain the true B^ lifetime by defin- 
ing a factor K, where ct = Kct* . We evaluate the K 
factor distribution, H(K), for £?+ events using Monte 
Carlo simulation. We are then able to express the 
distribution of ct* for B+ — > J/4>l + X as 



K Krt* 
F B (cf.tr) = / dKH{K)—0{ct*)exp{ ) ® G(o) 

CT CT 

(3) 

where cr is the average B+ lifetime, a represents the 
estimated error on the measurement of ct* for each 
event, and G(o~) is defined as 



G(a) = 



1 



2nsa 



(4) 



'Speaker, on behalf of the CDF Collaboration 



The measurement of the £?+ average lifetime is be 
carried out by minimizing —2Log(L), which is evalu- 
ated for the candidate J/ip + 1 events. L is the likeli- 
hood function for the ct* and a measured in candidate 
events and includes cr as a free parameter. 
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3. The B + cross section measurement 

c 

concept 

We measure 

a(B+)BF(B+ -+ J/^+u) = N(B+) 
a{B+)BF{B+ -» J/iPK+) N{B+) rel [ ' 

The basic strategy of the B+ cross section measure- 
ment is to reconstruct the number of £? + — > J /ip\x + v 
relative to the number of B + — > J /ipK + candidates, 
N(B+)/N(B + ), determine the relative detector and 
reconstruction efficiency, e re i — e(B + )/e(B+), and use 
these to determine the ratio of the final production 
cross section times branching ratio. 



4. The event selection 

The analysis presented here is based on the events 
recorded with a di-muon trigger that is dedicated to 
J/ip — > decays. Both analysis use the same 

J/ip — > n + [J.~ sample. The muon pair is reconstructed 
within a pseudo-rapidity range |r;| < 1.0. We select 
about 6.9xl0 6 J/ip candidates, measured with a mass 
resolution of approximately 12 MeV/c 2 . The di-muon 
invariant mass distribution is shown in Fig. 1. 
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Figure 1: The J/ip — > invariant mass distribution. 

Events within a mass range of ±50 MeV/c 2 around the 
central J/ip mass value were used for both the lifetime 
and cross section analysis. 

In addition to the di-muons from the J/ip decay, 
we require a third track that is matched to the same 
vertex as the J ftp- This third track can be from any 
of three samples of interest: £?+ — * J/ipl + ~K decays, 
B + — > J/tpK + decays, or just J/ip + track decays. 
The last sample represents sources of backgrounds to 
the B+ semileptonic decay. 



5. The B+ backgrounds overview 
5.1. Common for both analysis 

Both the lifetime and cross section analysis have 
some common sources of background: misidentified 
J/ip, misidentified third muons, and bb backgrounds. 
The misidentified J/ip background occurs when one 
of the muons is actually a mis-reconstructed hadron 
or muon from other sources that produce a mass con- 
sistent with that of the J/ip. The misidentified third 
muon background can arise from the following sources. 
The J/ip in the J/ip + track system is highly popu- 
lated by non-B+ sources. The third track associated 
with the J/ip could be a ir + or K + that can either 
decay-in-flight or punch-through the calorimeter and 
the steel absorber and produce the muon signature. 
The bb events represent cases when a J/ip is produced 
from one b jet and the third muon originates from the 
other b in same event. 

Figure 2 shows the misidentified muon rates for , 
K^- , and p(p) and the pseudo-proper decay length ct* 
for misidentified third muons. 




Hadron Pj (GeV/c) Pseudo-Proper Decay Length (|im) 



Figure 2: The misidentified muon rates from 7r , K , and 
p(p) as a function of hadron pr (left), and the pseudo- 
proper decay length ct* for misidentified third muons 
(right). 



5.2. Specific for each analysis 

There is an additional background for the cross sec- 
tion measurement. The selected Bj — > J/ip(i + X sam- 
ple contains contributions from other B^ decays with 
a tri-muon in the final state. For example, a B+ can 
decay into ip(2S)^ + u followed by ip(2S) -> J/ipX. 

The following backgrounds are specific to the life- 
time analysis: misidentified e^, residual conversions, 
and prompt J/ip. The misidentified e can arise from 
cases when a , K , or p from the J/ip + track 
system satisfies the likelihood function based on 
the calorimeter responses. The residual conversions 
are e from 7-conversion or tt° Dalitz decays. The 
prompt J/ip are additional J/tpl candidates where 
the J/ip originates from prompt non-i?+ sources. Fig- 
ure 3 illustrates the e + e~ veto efficiencies and the 
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pseudo-proper decay length ct* distribution for the 
J/^+Conversion e background sample. 
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Figure 3: The e + e~ veto efficiencies as a function of elec- 
tron pr (left), and the pseudo-proper decay length ct* dis- 
tribution for the J/^+Conversion e background sample 
(right). 



6. The B+ lifetime results 

6.1. Lifetime systematic uncertainties 

The systematic uncertainties in the -B+ lifetime 
measurement originate from uncertainties in our mod- 
els for background and signal events. Some of the 
largest systematic uncertainties are summarized be- 
low: 

• Resolution function - choice of model for detec- 
tor resolution: ±3.8 fim 

• Pythia model for bb background - relative con- 
tribution of QCD processes: ±2.4 fim 

• Vertex detector alignment - uncertainties in the 
positions of silicon detectors: ±2.0 fim 

• e + e~ veto efficiency - uncertainties related to 
modeling e + e~ veto efficiencies: ±1.5 /ira 

• B c spectrum - variations of the K factor distri- 
bution due to variations in the B c production 
spectrum: ±1.3 fim 

We add the individual uncertainties in quadrature 
to obtain a total uncertainty of ±5.5 fim. 

6.2. Lifetime results 

We fit the ct* distributions for signal candidates 
in the J /%pfi + and J/ipe + channels separately using 
likelihood functions based on our models for signal 
and background events. The fitted data is shown in 
Fig. 4 for B+ -> J/tp^+X and for B+ -> J/ipe+X 
decays. 

The B+ lifetime is found to be 179.1" 



L -27^( stat ) P- m 

for the J/ip^ + final state and 121.7+^3 (stat) fim for 
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Figure 4: The pseudo-proper decay length ct* distribu- 
tions for B^ — » J/ipfi + X decays (left) and for — > 
J/ifje + X decays (right) with their the background models 
superimposed. 



the J/ipe + decay mode, respectively. We performed 
the simultaneous fit of both samples and found an 
average B+ lifetime of 142.51^' ®(stat)±5.5(syst) fim. 
Figure 5 shows our average lifetime comparison 
with other measurements. 

CDF Run II Preliminary: ~1 fb" 1 
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Figure 5: World average of lifetime, which includes 
the CDF Run I B c lifetime, the most recent DO Run II re- 
sult, and the result presented in this paper. The lifetimes 
are weighted by the total variance of the individual mea- 
surements in the average. 



7. The B+ relative cross section 

In order to measure the -B+ relative cross section 
we need to find the numbers of B+ and B + , N(B+) 
and N(B + ), and determine the relative efficiency, e re i 
= e(B+)/e(B+). We select 229 (214) B+ candidates 
with the requirement pt(B^) > 4 (6) GeV/c, re- 
spectively. The number of -B+ signal events after 
backgrounds subtraction is presented in the follow- 



ing subsection. The number of B^ 



J/ipK + sig- 



nal events is found to be 2333 ± 55 (2299 ± 53) for 
Pt(B + ) > 4 (6) GeV/c, respectively. The combina- 
toric and B + — > J/ipir + contributions are subtracted. 
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7.1. The B+ backgrounds and excess 

The backgrounds and the resulting number of sig- 
nal events for the —>■ J /i]jfi + v decays are summa- 
rized in Table I. The background identified as "Dou- 



events are passed through the full detector and trigger 
simulation. The Monte Carlo simulation samples were 
processed in the same way as for the data. The effi- 
ciencies e B + and e B + for pt{B) > 4 (6) GeV/c, along 
with the relative efficiency, are presented in Table II. 



Table I Observed N(B+ -> J/ipfi+v) for the p T (J/ipfJ.) > 
4 GeV/c (6 GeV/c) threshold. 





p T {B) > 4 GeV/c 


p T (B) > 6 GeV/c 


N(B+) observed 


229±15.1(stat) 


214±14.6(stat) 


Misidentified J/ip 


21.5±3.3(stat) 


20.5±3.2(stat) 


Misid. third muon 


55.8±2.0(stat) 


53.6±1.9(stat) 


Doubly misid. 


-8.8±0.4(stat) 


-7.5±0.3(stat) 


bb background 


37.7±7.3(st+sys) 


35.4±7.0(st+sys) 


Other modes 


5.2±0.5(stat) 


4.8±0.4(stat) 


Total background 


111.4±8.3(stat) 


106.9±8.0(stat) 


signal 


117.6±17.2(stat) 


107.1±16.7(stat) 



bly misidentified" in Table I represents the subsam- 
ple of misidentified J/ if; and misidentified third muons 
that needs to be subtracted only once to avoid double 
counting. 

Figure 6 shows the invariant mass distribution of 
the — > J/i/j^ + X data events with the experimen- 
tal backgrounds and a Monte Carlo simulation of the 
signal sample superimposed (left), and with the back- 
ground subtracted (right). 



CDF Run II Preliminary: ~1 fb 1 CDF Run II Preliminary: ~1 fb 
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Figure 6: The invariant mass distribution of the B^~ — » 
J/ipfj> + X data events with the experimental backgrounds 
and a Monte Carlo simulation of the signal sample super- 
imposed (left), and with the background subtracted ap- 
plied (right). 



7.2. The relative efficiency e ret 

In order to determine the relative efficiency, we sim- 
ulate B+ -> J/iPK+, B+ -> J/ip^+v, and B*+ -> 
5+ 7 decays. As a description of the rj — pt spectrum 
for we use the most recent theoretical work by 
Chang et al. [3]. For the B + we used the spectrum 
from Ref. [7], which is found to be in good agreement 
with CDF measurements. All Monte Carlo simulation 



Table II Efficiencies for B+ and B + for p T (B) > 4 (6) 
GeV/c. 



Efficiency 


p T {B) > 4 GeV/c 


p T (B) > 6 GeV/c 


e Bt (%) 


0.0551 ±0.0010 


0.1232 ±0.0024 


e s+ (%) 


0.3231 ± 0.0022 


0.6005 ± 0.0042 




5.867 ± 0.068 (stat) 


4.873 ± 0.060 (stat) 



The pt spectra for data and Monte Carlo simulation 
are shown in Fig. 7. 



CDF Run II Preliminary: ~1 fb ' CDF Run II Preliminary: -1 fb" 1 




p T (JA|/+H*) (GeV/c) p T (JA|/+K*) (GeV/c) 

Figure 7: The comparison of the pt spectra of data versus 
Monte Carlo simulation for B^ — > J/ip(j, + i> (left), and for 
B + -> J/ipK + decays (right). 



7.3. Cross section systematic 
uncertainties 

We divide the systematic uncertainties into two cat- 
egories: uncertainties on the number of B+ signal 
events and uncertainties in the determination of the 
relative efficiency. 

7.3.1. B$ background uncertainty 

The systematic uncertainty considered in the de- 
termination of N(B^) arise from events which do not 
originate from decays. All backgrounds are as- 
signed a systematic uncertainty except the misiden- 
tified J/tp background. Since the estimation of this 
background is determined directly from the sidebands 
of the J/ip, we do not know of any source of systematic 
error that should be included. 

The largest source of uncertainty in the misidenti- 
fied third muon calculation is due to the poor knowl- 
edge of the proton fraction in the J/ip + track sample. 
The particle identification method, dE/dx, does not 
allow us to separate protons from kaons in our kine- 
matic region. Consequently, we measure the proton 
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fraction in other momentum ranges using the time- 
of-flight (TOF) particle identification combined with 
dE/dx information and then extrapolate the fraction 
to our momentum range according to measured trends 
in Monte Carlo simulation. 

The other N(B+) systematic uncertainty arises 
from a poor knowledge of non-exclusive £>+ — ► 
J/ipp, + X branching ratios and is estimated by vary- 
ing the branching ratios of eleven B~£ decay modes 
that may contribute to the sample of tri-muon events. 
In order to assign a systematic uncertainty we double 
and halve the branching ratios of the non-exclusive 
decays with respect to the rate of £>+ — > J/ipp, + v. 
We choose the larger variation in cither direction for 
the systematic uncertainty. 

Table III summarizes all of the background sys- 
tematic uncertainties assigned. 

Table III Systematic uncertainties on the number of 
— > J /ipu + v events for different pr(B) thresholds. The 
total uncertainty is calculated by adding all the individual 
uncertainties in quadrature. 



N(B+) uncertainties 


p T (B) > 4 GeV/c 


p T (B) > 6 GeV/c 


Misid. third muon 


±5.7 (sys) 


±5.5 (sys) 


Doubly misid. 


±0.9 (sys) 


±0.8 (sys) 


Other Be decays 


+ - 6 2 °s (^) 


tS:S (sys) 


Total 


±U (sys) 


tli (sys) 



• the combined £?+ + B* + spectrum and a pure 
B + spectrum. 

The systematic uncertainty due to knowledge of the 
B + pt spectrum is estimated by varying the Monte 
Carlo simulated spectrum below 10 GeV/c to bring 
it into agreement with the data (see the right plot 
in Fig. 7. The difference between the nominal and 
recalculated relative efficiency is assigned as the un- 
certainty due to B + pt spectrum. 

Another source of systematic uncertainty that we 
consider is the different XFT efficiencies of kaons and 
muons that exist in the data and are not modeled in 
the simulation. 

The total e re i systematic uncertainty is summarized 
in Table IV. 

Table IV Systematic uncertainty assigned to e re i for dif- 
ferent p T (B) thresholds. The numbers 0.720 and 0.298 
in the "Total" represent the uncertainties due to of 
spectrum. 



e re i uncertainties 


p T (B) > 4 GeV/c 


p T (B) > 6 GeV/c 


_B+ lifetime 


+0.393 
-0.223 


+0.354 
-0.160 


B^ spectrum 


± 0.720 


± 0.298 


B + spectrum 


± 0.340 


± 0.161 


XFT systcmatics 


± 0.192 


± 0.160 


Total 


±S;SS±0.720 


1° ;? 7 °±0.298 



7.3.2. Relative efficiency systematic uncertainty 

We consider the systematic uncertainty in the pre- 
diction of the relative efficiency due to the measured 
statistical uncertainty of the -B+ lifetime, knowledge 
of the production spectra for and B + , and differ- 
ences between K and fi triggering rates at the first 
level of the CDF trigger system, the extremely fast 
trigger, XFT. 

The relative efficiency systematic uncertainty due 
to B+ lifetime uncertainty is estimated by varying the 
lifetime by ±14 /im of its default value. This variation 
represents approximately one standard deviation in 
the average lifetime result. 

The relative efficiency systematic uncertainty due 
to the knowledge of £>+ pt spectrum is fully based 
on the theoretically predicted pr spectra from work 
in Ref. [3]. We consider the variations between: 

• doubling the qq contribution relative to the nom- 
inal approach; 

• the pure gluon fusion model, called the fixed fla- 
vor number ( FFN) model, and the more com- 
plete gg + gb + gc model, known as the general- 
mass variable-flavor-number (GMVFN) model; 



8. The B+ relative cross section results 

We have performed a measurement of the relative 
production cross section of £?+ — * J/ip^iy in inclu- 
sive J/ip data with an integrated luminosity of 1 fb _1 . 
We have identified a sample of 229 (214) events with 
an estimated background from all sources of 111±8 
(107±8) events for p T (J/iPp+) > 4 (6) GeV/c, re- 
spectively. The final numbers used in the cross section 
measurement, including systematic uncertainties, are 
given in Table V. 

We give the result for the ratio 
ts^CZi:] with M B) > 4 GeV/c 
thresholds as 

0.295 ±0.040 (stat.)+°;°2g (syst.) ± 0.036 (p T spec) 

and for pr(B) > 6 GeV/c as 

0.227 ±0.033 (stat.)+o:oi7 (syst.) ±0.014 (p T spec). 

Of the two results, the measurement with the 
Pt(B) > 6 GeV/c threshold has the lower systematic 
error. Below 6 GeV/c there is uncertainty in the B + 
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Table V Final numbers used in the calculation of the rel- 
ative Be — > J/ipu + v production cross section times the 
branching ratio for two different pr{B) thresholds. 



Final values 


p T (B) > 4 GeV/c 




If7.6±f7.2 (stat) tgj(sys) 


N(B+) 


2333± 55 (stat) 




5.867±0.068 (stat) to^sys) 
±0.720 (B+ spectrum) 


Final values 


p T (B) > 6 GeV/c 


N(B+) 


107.2±16.7 (stat) ±™ (sys) 


N(B+) 


2299±53 (stat) 




4.872±0.060 (stat) I^'Its ( s y s ) 
±0.298 (B+ spectrum) 



efficiency that appears to introduce a significant sys- 
tematic discrepancy between the simulated spectrum 
and the spectrum as determined from the data. 

Using theoretical assumptions and independent 
measurements, we are then able to calculate the to- 
tal cross section. Using the measured quantities 
BR(B+ -> J/ipK+) = (1.007 ±0.035) x 10~ 3 [8] and 
a(B+) = 2.78 ± 0.24 ^b for p T (B+) > 6 GeV/c [4], 
we calculate 

a(B+) * BR(B+ -» J/ipLi+is) = 0.64 ± 0.20 nb 

ior Pt(B+) > 6 GeV/c. Assuming that the branching 
ratio BR(B+ -> J/^ii+v) = 2.07 x 10~ 2 [9], we find 
the total cross section to be 

a{B+) = 31 ± 10 nb 



9. Conclusions 



We have performed measurements of the life- 
time and production properties based on semileptonic 
£>+ — ► J/ip ± l + + X decays using data from pp col- 
lisions collected with the CDF II detector correspond- 
ing to an integrated luminosity of 1 fb _1 at y/s=l.96 
TeV. 
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